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Abstract 

We investigate the dynamics of a nonequilibrium interface between coexisting 
phases in a system described by a Cahn-Hilliard equation with an additional driving 
term. By means of a matched asymptotic expansion we derive equations for the 
interface motion. A linear stability analysis of these equations results in a condition 
for the stability of a flat interface. We find that the stability properties of a flat 
interface depend on the structure of the driving term in the original equation. 
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1 Introduction 



Off-equilibrium systems composed of regions with different phases can exhibit a variety 
of patterns, such as fingers aligned in a certain direction. A familiar example is spinodal 
decomposition of binary mixtures (for a review see ||, @). When quenched below its criti- 
cal temperature an initially homogeneous system develops domains of the new equilibrium 
phases. After some time these domains are well defined, i.e. they are separated by sharp 
interfaces. The system further approaches equilibrium by means of the motion of these in- 
terfaces, the initial pattern evolves in the course of time. Another example is the dynamics 



of driven diffusive systems (for a review see Ref. [0). Interacting particles undergoing 
biased diffusion tend to form clusters under certain conditions. Like for the first example 
there are well defined interfaces, which undergo some dynamics. In all of these phase segre- 
gating systems patterns are selected according to their stability. For example the instability 
of a flat interface causes the growth of fingers out of an initially flat interface. For this 
reason one would like to understand the stability properties of nonequilibrium interfaces. 

On a course grained level these phase segregating systems are described by means of 
an order parameter, such as a particle density. Typically, the dynamics of this order 
parameter can be modeled by means of some nonlinear differential equation. One of the 
simplest such equation is the Cahn-Hilliard equation (CH equation), which was introduced 
in the context of binary mixtures [|J. The properties of this equation are rather well 
understood. In particular, the dynamics of the domain boundaries obey a set of linear 
differential equations, which have been derived by Pego || using of a matched asymptotic 
expansion. An instability of the Mullins-Sekerka type has been observed. 

Certain modification of the CH-equation were introduced to take into account the effect 
of external fields PH-fllOl, Hl2 | (see Refs. [fTl|, |13| , |14| for modified CH equations in the 



context of driven diffusive systems). One of the questions one would like to answer is the 
following: How does an external field change the stability properties of a flat interface? 
This question was addressed by Yeung et. al. [IDj| for a certain type of driving term. This 



driving term has the form of an additional current which is proportional to the field and 
a conductivity. The conductivity was assumed to be a concave second order function of 



the order parameter with the same symmetry as the chemical potential. Yeung et. al. J10 
found that the field modifies the interface dynamics. An interface perpendicular to the 
external field is stable for one direction of the field and unstable for the other direction. 

In this paper we generalize the work of Yeung et. al. |K|. We assume a more general 
structure of the driving term, i.e. a more general conductivity. Our approach is the same 
as the one from Ref. |T(| (see also Ref. 0). Assuming that the driving field is sufficiently 
small we perform a matched asymptotic expansion of a modified Cahn-Hilliard equation 
with an arbitrary conductivity. The resulting set of linear equations describes the motion 
of an interface between coexisting phases. With the help of these equation we analyze the 
stability of a flat interface. We find a general stability criterion, which takes a very simple 
form if we assume that the conductivity is a third order polynomial in the order parameter. 
For such a conductivity we show that an interface perpendicular to the external field can 
be unstable for either direction of the field, and that an interface parallel to the external 
field is always stable. 

The paper is organized as follows. In Section 2 we define the model. A matched 
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asymptotic expansion is done in Section 3. Equations for the motion of an interface are 
derived. In Section 4 we perform a linear stability analysis for a flat interface resulting in 
a condition for the stability of such an interface. This condition is subsequently shown to 
take a simple form for the case of a conductivity that is a third order polynomial in the 
order parameter. Section 5 summarizes the results. 



2 The model 

We consider a two-dimensional system described by a conserved order parameter iz(R, T) 
where R denotes a position in space and T the time. The dynamic is assumed to be defined 



by the following modified Cahn-Hilliard equation @-||, [fLOl , |[T2| , jnj: 

% ~ — w 

J = -V/i + E a{u) (2) 
[i = \u{u - - M m ) - £ 2 V 2 w . (3) 

In the above equations, J is the current of the order parameter, fi is the chemical potential 
of the free system, E is an external field (e.g. an electric field), and cr(u) is the conductivity. 
One may think of the chemical potential as related by means of fi = 5F/5u to a free energy 
functional F[u] = J dR { f(u(H)) + £ 2 |Vm(R)| 2 /2 } , where f(u) is a bulk free energy 
density the derivative of which is f'{u) = \u{u — ^)(w — u m ) . The current has two 
additive contributions: One that minimizes the free energy and another one defining the 
transport caused by the field. The first one is specified by the structure of the chemical 
potential. Assuming A > and u m > Eq. @ defines a chemical potential such that the 
free energy has minima in regions with u = or with u = u m > 0. The term proportional 
to £ 2 results in an increase of energy whenever there is a gradient in u, i.e., £ 2 stands for 
a surface tension. The current caused by the field is E cr(u) , it depends on the function 
a(u) . The standard conductivity used in the literature ||-[f§, [[l(J, 13; @| is a concave 
second order function of u , which has the same symmetry as the chemical potential, i.e. is 
symmetric around u m /2. Here, we do not restrict to this type of function, but rather want 
to see what kind of behavior can be observed for a more general conductivity. 

Equations of type ([]])- (|3|) are usually motivated on a purely phenomenological (see e.g. 
0, H, [10|, [12], |13||) or mean field (see e.g. 0) level. The standard type of argument results 
in a mobility as a factor to the total current, i.e. a current of the type a{u)[— V/i + E] . 
Assuming furthermore that the variation of cr(u) is small enough cr(w)V/i is then replaced 
by V/i . We follow the same logic here. It is also worthwhile mentioning that a rigorous 
derivation of a more general type of equation has been done for diffusive system with 
long-range interaction and Kawasaki dynamics [ |TT| . 

Since it is more convenient to work with a dimensionless equation we rescale Eqs. 
(S)"©- Simultaneously we add a constant term to the conductivity and apply a Galilean 
transformation such that the new current vanishes at both minima of the chemical potential. 
The adding of a constant does not change the equation for u(R, T) , since only the gradient 
of a enters into this equation. The Galilean transformation is nothing but a change of the 



3 



reference frame. We define a length Lq and a time Tq obeying 



1 = L 2 T Xu 2 m (4) 

as well velocity 

V g = |E| u m l a(u m ) . (5) 
Now we define dimensionless quantities as: 

(6) 
(7) 
(8) 





r = 


L^R-L^VgT e 




t = 






P = 




where 




E 

e = i — : 
|E| 



(9) 

is the unit vector in the direction of the field. Inserting the above definitions into Eqs. 
(HI)"© yields the following dimensionless equation 

| - -v.j do) 

j = _ V z/ + e/?(p) (11) 
v = p(p-l)(p-l)-e 2 V 2 p (12) 



with 
and 



e 2 = LZT\- L u^ (13) 



(3(p) = e^mX~ 3 / 2 u m 4 { a{pu m ) - a(0) - p[a(u m ) - a(0)] } . (14) 

In Eqs. (PH|)-(rr^) the operator V stands for derivative with respect to the new coordinate 
r . The conductivity (3{p) in the new frame has the property 

(3(1) = 0(0) = (15) 

i.e., it vanishes at both minima of v . However it is not necessarily positive in the whole 
interval (0, 1) . In fact, later on we will identify a regime where f3(p) < for some part of 
(0, 1) . For that reason it would be hard to give a direct physical meaning to Eqs. (PH)-(PD 
in terms of the rescaled variables. 

The velocity V g given by Eq. (|5|) is positive by definition. Consequently the new frame 
of reference moves in the direction of the field with respect to the original one. 



4 



3 Dynamics of an interface 



The model denned in the previous section describes phase segregation, i.e., the order pa- 
rameter p(r, t) evolves such that there are regions in space where p(r, t) approximately 
takes the values for which the chemical potential is minimal. Since there are two minima 
of the chemical potential, there are two types of regions, which are separated by transition 
layers. The transition layer between the different regions, which is seen as a line on a large 
enough scale, evolves during time. We are going to give an analytical description of its 
motion for the case where 

e < 1 (16) 
and for < p < 1 : (3{p) = O(e ) and (3'{p) = 0(e°) (17) 

in the rescaled equations (|l(])-(0)- m terms of the original parameters these conditions 
read 

^A-^- 1 « 1 (18) 
and |E|f \- 3/2 u m 4 a max < 1 (19) 

where a max denotes the largest values a{u) and its derivative cr'{u) take in [0, 1] . The first 
of these conditions can be fulfilled for any choice of the parameters simply by choosing a 
large enough length scale Lq . However, the second condition poses a restriction on the 
parameters, it can be understood as a small field condition. 

Starting with the rescaled equations (|T0|)-(|T2|) we apply a matched asymptotic expan- 
sion. The basic idea is to expand the order parameter and its evolution equation in powers 
of the small parameter e. We closely follow the line of Refs. M and pi. Since we want 



to describe the motion of a single interface, we consider an initial configuration with two 
semi-infinite regions Q + and Q~ and assume that the order parameter is equal to 1 (0) in 
Q + up to corrections of order e . The values 1 and are the minima of the chemical 
potential. The regions Q + and Q~ are separated by a transition layer. We assume this 
layer to have a width of order e like it in the field free case ||, pi. The characteristic 
time of its motion is 

t = e t = e Tq 1 T . (20) 

In order to have a clear definition on any scale we define the interface T as the set of points 
where p = | . (Our results are not sensitive to the choice of the value 1/2 .) 

We are going to expand our equations on the time scale r first in the bulk of Q + and 
Q~ and then in the transition layer, i.e. near the interface T . 

3.1 Equations far from the interface 

Far from the interface T we define the new field 

p ± (r,r)=p(r,e- 1 t) (21) 

and similarly z^ ± (r, r) and /3 ± (r, r) where the superscript denotes the region we are consid- 
ering. Eqs. (|ToD and (0) result in the following equations on the time scale r 

ed T p ± (r,r) = VV(r, r) - e ■ V/^r, r) (22) 
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where v ± and (3^ are defined by Eqs. ([12]) and (|14D , respectively. Next we expand p ± (r, r) 
in powers of e as 

p ± (r,r) = p^ + epf(r,r) + 0(e 2 ) (23) 

where 

p+ = 1 and Pq = (24) 

according to our previous assumptions. Inserting above expansion into Eq. (|22|) , (|i~2|) and 
([T4]) results in an equation for pf(r,r) where the terms of zeroth order in e cancel each 
other. Comparing terms of first order in e yields the following relation 

[ V 2 -2B ± e- V] pf(r,r) = (25) 

with 

B+ = I3'{1) and B" = P'{0) . (26) 



3.2 Equations near the interface 

We are going to expand Eqs. (|I~0|)-(|T2"|) in the transition layer where the chemical potential 
is not close to its minima. Let us denote by r r a point there, i.e. a point with a distance of 
order e from T, and introduce the quantity $(rr,r) as the signed distance between point 
rr and the interface T at the time t = e _1 r. The sign is chosen such that $ > in Q + and 
$ < in Q~ . Furthermore we define 

m(r r ,r) = V$(r r ,r) (27) 
K(rr,r) = V 2 $(r r ,r) (28) 
V(tt,t) = 9 r $(r r ,r) . (29) 

As can be seen from these definitions, if r-p lies on T the vector m is the unit normal of 
the interface at r-p, k is the curvature there and V is the normal velocity of T on the time 
scale t. Signs are such that the vector m points into Q + , the curvature k is positive if the 
center of curvature lies in Q~ and V is positive if the interface moves towards Q~ 

Since the transition layer has a width of order e it is convenient to introduce a variable 

z = e- 1 $(r r ,r) (30) 

and for an arbitrary field /(rr, t) a field f(z, r r , r) by: 

/(r r ,t) = /(e- 1 $(rr),r r ,et) (31) 

The derivatives occurring in our equations are transformed as follows: 

V/ = [ V r + e-'mdz ] f (32) 
V 2 / = [ V 2 . + e-\d z + e~ 2 d 2 z ] / (33) 
dj = [e d T + V d z ] f (34) 

where V r is the gradient acting on rr only, and arguments were dropped for convenience. 
Assuming furthermore that all relevant fields near T depend only on their relative position 
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with respect to T, i.e., that f(z,rr + 5 m(rr, t),t) = /(z,rr,r) for any small number S, 
we find 

V$(r r ,r)- V r /(z,r r ,r) = . (35) 
With the above relations we are ready to transform and expand in powers of e as: 

P = P = Po + e pi + e 2 p 2 + 0(e 3 ) (36) 
v = v = v + ev 1 + e 2 v 2 + 0(e 3 ) (37) 
P = = P + e Px + e 2 p 2 + 0(e 3 ) (38) 

where the arguments, which we have dropped for convenience, have to be taken according 
to Eq. (|31f) . The expansions of v and a are related to the one of p by means of Eqs. dl2|) 
and (III). These relations can be casted in an explicit form using Eqs. (P^j), (§3|), (0) anc ^ 
(|35|). The following ones we are going to use later: 



^0 


= Po(Po-l)(Po- 


~\) - d 2 J 




(39) 


V\ 


= [3p (po - 1) + 


i] pi - Kd z p Q 


- fifo 


(40) 




= [3po(Po - 1) + 


1. 3, 
^ + 2( 2 Po 


- l )Pi - V r 2 p - Kd z p! - <9 2 p~ 2 


(41) 


Po 


= Pih) 






(42) 


Pi 


= PiP'(Po) ■ 






(43) 



The last two of these relations hold under the condition fllTD only. Inserting now expressions 



(]32|) , (p3D and (gj) as well as expansions (p6|) , ( pTj ) and (pq) into our evolution equation 
(10), (flT|) and comparing terms of order e~ 2 ,e _1 , e°, respectively, results in the following 
set of equations: 

= d% (44) 
= k d z u + d 2 ui — e ■ m d z p (45) 
Vd z p = V 2 z> + K d z v x + d 2 z v 2 - e ■ V r /3 - e ■ m <9 2 /3i . (46) 

The dynamics of the order parameter near the interface is described by the above equations 
combined with the following boundary conditions: The solution of above equations should 
match p ± (r r ,r) in fl ± , i.e., outside the transition layer. We demand 

lim p(z,Tr,r) = p ± (r T ,r) . (47) 

z — >±oo 

Using lim^-i-oo p(z, r r , r) = lim^-too p(r r + e z m, e _1 r) and expanding above condition 
in powers of e yields 



lim p (z,r r ,r) = p% (48) 

Z — »±oo 

lim pi(«,rr,r) = pf(r r ,r) (49) 

Z — »±oo 

lim p 2 (z,r r ,r) = pf(r r ,r) + lim [ z m • Vpf{r r , r) ] . (50) 

Z— >±00 2^±00 



Here we used Eqs. (gTj) , (|D and (||). 

One can now solve Eqs. (13), fl3TS|) with the boundary condition (^). The result is 



p (*,r r ,r) = [1+ e'^r 1 (51) 
u (z,r r ,r) = . (52) 

Let us now define the following integrals for later use: 

JdMz)?dz = ^ (53) 
°° [p (z)] n (3(p (z)) dz = V2 f 1 ^— (3(p) dp . (54) 

-oo JO 1 — p 

The first integral was computed and the second one was simplified using Eq. ([[I]). Next we 
insert 9 = into Eq. (f45|) and integrate over z . We fix the integration constants by means 
of the limits z — *> ±oo. In order to compute these limits we expand lim^-toc, 9(z, rp, r) = 
lim^^-too z/(rr + e z m, e"V) in powers of e and use the fact that /5(0) = f3{l) = . We 
obtain: 

d z v x = e • m & . (55) 

Integrating once more over z and taking again the limits z — * ±oo yields a relation between 
Pi and />f . Another relation is obtained by multiplying Eq. (f45"D by <9 z p and integrating 
over z from — oo to +oo. Combining those two yields 

p+(r r ,r) = - 2S , «(r r ,r) + 2I X e • m(r r , r) (56) 
pr(r r ,r) = - 25 , K(r r ,r) + 2(/i-/ ) e.m(rr,r) (57) 



where Jo and Ji are the integrals defined by Eq. fl54|) . 

In a last step we integrate Eq. ( |46| ) and take the limits z — > ±oo. As a result we obtain 
an expression for the interface velocity 

V(t t ,t) = ^m(r r ,r) [ Vp^(r, r) - Vp^(r, r) _ 

+ [ Jo + 2(£? + - ] e ■ m(r r , r) /<r r , r) 

- 2 [ (5+ - B")/! + B" J ] [e ■ m(r r , r)] 2 (5* 



where B± is defined by Eq. (||), S 1 by Eq. (H) and J , J x by Eq. (H). 

Above equations provide a macroscopic description of the interface dynamics. Macro- 
scopic means that the space resolution is of order e° so that the transition layer can be 
identified with interface T, i.e., any point rp can be considered as lying on T. Suppose 
that at a time t = e~V there is an interface T with a unit normal m(rr,r) and a cur- 
vature n(r-p, t) at its points rp. Away from the interface the order parameter is given by 
p ± (r r ,r) = Po + e pf( r r, T ) where pg = 1 , pg = and pf(v r ,r) is given by the solution 
of Eq. (^) with boundary conditions on the interface defined by Eqs. (0), (|)7]) i n terms 
of m(rr,T) and /t(rr,r). This solution in its turn determines the normal velocity of each 
interface-point by means of Eq. (|58|) . 
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4 Linear Stability Analysis 

As derived in the previous section, on the scale r = e _1 t = e _1 Tq 1 T the motion of 
an interface at r = rr separating the two regions Q + (with a particle density p + (r,r) = 
1 + epf(r, t) + 0(e 2 )) and Q~ (with p~(r, r) = + epi"(r, r) + 0(e 2 ) ) is given the normal 
interfacial velocity V of the interface 

V(r r ,r) = im(r r ,r) [ Vp+(r, r) - Vp^(r, r) ]|^ 

+ [ Jo + 2(5+ - £T)S ] e • m(r r , r) «(r r , r) 

- 2 [ (£+ - S-) A + B~ J ] [e • m(r r , r)] 2 (59) 

where m is the local normal to the interface, k = V • m is the local curvature, and B , S, Iq 
and Ji are constants that depend exclusively on the form of the conductivity a and are 
given by Eqs. (^),(Q) and fl5I|). The velocity depends on pf(r, r) which satisfy the linear 
partial differential equation 

[ V 2 -2B ± e- V] pf(r,r) = (60) 

subjected to the boundary conditions 

p+(r r ,r) = - 2S k{t t ,t) + 2h e-m(r r ,r) (61) 
pr(r r ,r) = - 2S k(t t ,t) + 2(h - 1 ) e-m(r r ,r) (62) 

on the points rp of the interface. 

We are interested in the stability of a flat interface against small perturbations in its 
profile. Let us consider an interface T of the form 



Vr 



(x,t) = he i{kx - WT) + V t ; kh « 1 (63) 



(see Fig. 1) separating two semi-infinite regions Q + (where y > yr{x,r)) and Q~ (where 
y < yr{x, r)). The direction e of the external electric field is arbitrary for the moment. To 
first order in kh, we have 



-ikhe i( - kx -^ and m y = 1 (64) 



for the x and y components of m respectively. The local curvature k of the interface is 
given by 

K = k 2 he i( - kx ~^ (65) 
while the normal velocity V takes the form 

V = -dy s /dT = -V + iuhe i( * kx ~ wr) . (66) 

Eq. (|60|) has a solution pf (r, r) of the form 

pt(v,r) = ^(C) + kh e^ kx -^g ± (0 + O(kh) 2 (67) 
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with £ = y — he 1<ykx The boundary conditions Eq. (|6l| ) and (|62"D are satisfied for 

/ + (C) = 2he y -C+ (l-e 2B+ ^) (68) 

g +(() = -[2S k + 2ie x h]e- x+ <-2C + ^ ( e ~ A+ < - e 2B+e ^) (69) 

/-(C) = 2{h-I )e v -C- {l-e 2B ~^) (70) 

<T(C) = -[2S k + 2ie x (h - I )]e- X ~ ( - 2C~^ (e~^ - e 2B ~^) (71) 

where e x and e y are the x and y components of the unit vector e pointing in the direction 
of the external field E , and A are solutions of the following quadratic equation: 

(A ± ) 2 + 2B ± e y X ± - k(k + 2ie x B ± ) = (72) 

subjected to the constraint Re(X + ) > (Re(\~) < 0). The new constants C + and C~ are 
determined by the boundary conditions at z = ±oo. We will discuss the effects of (7 ± in 
the following section. 

Inserting pf (r, r) given by Eqs. ( |5"7| ) - ( |7T| ) in the equation for the normal velocity Eq 
(|59|), we can compare with Eq. (|66"D to obtain an expression for the constant velocity Vq 

V = B + (2he y - C + )e y - B~ [ 2{I X - h)e y - C~\e y (73) 
= [B + p+(z = +oo) -B~p{{z = -oo)] e y (74) 

and for ui 

iu = hj-(g + (z)-g-(z))\ z=0 + k 2 [I + 2S(B + -B-)}e y 

+ Ak [I B~ + h(B + - B~)] ie x e y (75) 

Assuming that are real, and writing X ± = Xf + iXf and u> — u r + iu>i, we get 

Ui = -k 2 [I Q + 2S(B + - B~)] e y - k 2 S(X+ - X~) 
+k [hXf - (h - I )Xr] e x 

-C + B + (X+ + 2B + e y )e y + C~B~{X; + 2B~e y )e y (76) 



A flat interface is unstable against small perturbation of the form ( |63"D if the external 
field, the conductivity and the boundary conditions at infinity are such that u>i > 0. Eq. 
is the main result of this paper. 



4.1 The case of a third order conductivity 

In order to further study the properties of Eq. ( [76]) we assume in this section that the 
conductivity o{u) can be written as a polynomial of the form 

cr (it) = cr + V\U + a 2 u 2 + cr 3 u 3 . (77) 

Previous studies pil of the stability of a flat interface in a driven Cahn-Hilliard system of 
the form Eqs (]!])- (H) have limited themselves to the case where a{u) = u{u m — u). While 
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we can reproduce their results, we will observe interesting new effects considering the more 
general form ( [77|) of a(u) . 

The constants in Eq. (|76|) can be expressed explicitly in terms of the system parameters. 
Using Eqs. fl2fj), (^), Q5§ and (@) we find: 

5+ = A(a 2 + 2a 3 M m ) (78) 
fi- = -A (a 2 + a 3 u m ) (79) 

Jo = -V2A(a 2 + ^a 3 u m ) (80) 

1 5 

A = -V2A (-a 2 + -a 3 u m ) (81) 

2 o 

S = ^ (82) 
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with 



A = e - 1 |E|eA- 3/2 M - 2 > . (83) 



All of the above constants are independent of do and o\. The reason for this is that the 
equations for the order parameter do not depend on o"o, and that a change of o\ amounts 
to a change of the reference frame only 

4.1.1 External field parallel to interface 

In this case the direction of the field is such that e x = ±1, e y = 0. To lowest order in k we 
get, using Xf « ±k 1 l 2 (\B ± \) 1 / 2 and Xf = 

m = k" 2 [h^= + (h - Io)-^=] ■ (84) 



Using the explicit expressions (|Z8|) - (|83|) one can show that u { < 0, i.e. that the interface 
is stable, for any choice of {cr , <Ji, cr 2 , 03}. This result does not depend on the sign of e x , 
nor does is depend on the actual boundary conditions away from the interface through the 
constants C ± . 

4.1.2 External field perpendicular to interface 

In this case we have e x = and e y = +1 if the external field points into the high density 
region Q + or e x = and e y = —1 if it points towards the low density region Q~. We also 
have Xf = from Eq. (0). 

Depending on the signs of B + e y and B~e y we get different conditions on the possible 
values of the constants C ± . From Eqs ( ]S5| ) we see that if B + e y > (B~e y < 0) we must 
take C + = (C~ = 0) in order for pf (pj) to be finite as z — >• 00 (z — >• —00). This means 
that the system does not support a single interface in this case unless there is some current 



of particles at infinity, (see also the discussion of the one-dimensional model in Ref. ||12| ) 
If B + e y < (B~e y > 0) then there are no such constraints on C + (C~) and we chose 

the special cases C + = 2I\e y (C~ = 2(I\ — h)e y ) which makes pf (p7) vanish at infinity. 
We are left with four cases to study for uf. 
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CASE A: B + e y > and B~e v < 0. If we assume that k « B ± , then Xf « 



and we get, to lowest order, 



-k 2 e y [Jo + 2S(B + - B-)} 



(85) 



CASE B: B + e y < and B~ e y > 0. To lowest order we write A^ = -2B ± e y (1 + 

4(B±)2 _ 16(B±)*)' alld § et 



/c 67, 



A(B+) 2 (B-) 2 



[h{B-f - (h - I ){B + f + 2S(B+ - B-)B + B-} (86) 



CASE C: B+e y < and B~e y < 0. With A+ = -2B+e y (l + j^+p), A; 
we get 



/C €.y 

' 2B- ' 



(87) 



CASE D: -B + e ?/ > and B e y > 0. This is essentially the reverse of the previous 



case. Aj!~ — , A r 



4(b-)2 )> an d get 
-k 2 e y (h + 2S B+) 



(88) 



Inserting the explicit expressions (|78|)-(|8"3|) for B ± , I , I\ and S into above equations 
for Ui one comes to the following conclusion: In all four cases the interface is unstable if 



(2cr 2 + 3a 3 u m )e y > 



(89) 



and stable otherwise. Condition (|89D is the main result of this section. 

It is natural to assume the conductivity to be positive in the whole interval [0,w m ] . 
Assuming furthermore that cr(u) is symmetric around u m like the chemical potential and 
vanishes for u = , we are lead to the standard expression a(u) = u(u m — u) There 
o"2 < and C3 = , the interface is unstable if e y — — 1, i.e. the case where the external field 



points away from the high density region Q + [ 10[ . Although the assumptions of symmetry 
and vanishing <r(0) are natural for certain Ising-like ||-||, \TU\, [JEJ, |T3[ or particle-hopping 



models | llf , one can think of more general situations (e.g. the microscopic model in [p 
Then, the parameters {a , G\, a 2 , 03} can be such that (2cr 2 + 3o- 3 u m ) > . In this case 
the interface is unstable for e y = +1 , i.e. if the external field points into the high density 
region 



5 Summary 

We have applied a matched asymptotic expansion to a system of equations describing the 
dynamics of phase segregation in the presence of an external field. The influence of the 
field on a region of local density u is given by a conductivity a{u). We derived equations for 
the dynamics of an interface separating two regions in different phases (see Eqs. (|25|) , ([56]), 
( |5T| ) and ( |58"D ) and studied the stability of a flat interface against small perturbations. We 
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found a general condition (see Eq. (|76|) ) for the stability of such an interface and discussed 
in more detail the case where u{u) can be written as a third order polynomial in « . In 
this case the interface is always stable if the field is parallel to it. However, if the field is 
perpendicular to the interface, the interface can be either stable or unstable depending on 
the conductivity and on the direction of the field (see Eq. (|^) ). 



Acknowledgements 

We thank B. Schmittmann and R. K. P. Zia for stimulating discussions. S. S. gratefully 
acknowledges financial support by the Deutsche Forschungsgemeinschaft. C. A. L. has been 
supported by the National Science Foundation through DMR-9419393. 

References 

[1] J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958) 

[2] J. D. Gunton, M. San Miguel and P. S. Shani, in Phase Transitions and Critical 
Phenomena, Vol 8, C. Domb and J. L. Lebowitz, eds., Academic Press, New York 
1983 

[3] J. S. Langer, An Introduction to the Kinetics of first- order Phase Transitions, in Solids 
far from Equilibrium, C. Godreche, ed., Cambridge University Press (1991) 

[4] K. Kitahara, Y. Oono and D. Jasnow, Mod. Phys. Lett. B 2, 765 (1988) 

[5] J. S. Langer, N. Bar-On and H. D. Miller, Phys. Rev. A 11, 1417 (1975) 

[6] S. Puri, N. Parekh and S. Dattagupata, J. Stat. Phys. 75, 839 (1994) 

[7] A. M. Lacasta, A. Hernandez- Machado and J. M. Sancho , Phys. Rev. B 48, 9418 
(1993) 

[8] A. J. Bray, Adv. Phys. 43, 357 (1994) 

[9] R. L. Pego, Proc. R. Soc. Lond. A 422, 261 (1989) 

[10] C. Yeung, J. L. Mozos, A. Hernanez-Machado and D. Jasnow, J. Stat. Phys. 70, 1149 
(1993) 

[11] G. Giacomin and J. L. Lebowitz, Phys. Rev. Lett. 76, 1094 (1996) 
[12] C. L. Emmot and A. J. Bray, Phys. Rev. E 54, 4568 (1996) 

[13] F. J. Alexander, C. A. Laberge, J. L. Lebowitz and R. K. P. Zia J. Stat. Phys. 82, 
1133 (1996) 

[14] B. Schmittmann and R. P. K. Zia, Statistical Mechanics of Driven Diffusive Systems, 
in Phase Transitions and Critical Phenomena, vol. 17, ed. by. Domb and Lebowitz, 
Academic Press, U.K. (1995) 



13 



[15] L. K. Wickham and J. P. Sethna, Phys. Rev. B 51, 15017 (1995) 



Caption to the figure 

Fig.l: Schematic picture of the approximately flat interface V defined by Eq. flB3|). The 
interface separates the semi-infinite regions Q + and fl~ , in which p sa 1 and p ~ , 
respectively. In the figure m denotes the unit normal vector on T and h the amplitude of 
the perturbation around a flat interface. The unit normal vector of the external field is 
denoted by e. 
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